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Abstract 

We consider the complete probability distribution P({T„}) of the trans- 
mission eigenvalues Ti , T2 , . . . T/v of a disordered quasi-one-dimensional con- 
ductor (length L much greater than width W and mean free path /). The 
Fokker-Planck equation which describes the evolution of P with increasing 
L is mapped onto a Schrodinger equation by a Sutherland-type transforma- 
tion. In the absence of time-reversal symmetry (e.g. because of a magnetic 
field), the mapping is onto a free-fermion problem, which we solve exactly. 
The resulting distribution is compared with the predictions of random-matrix 
theory (RMT) in the metallic regime (L <C Nl) and in the insulating regime 
{L ^ Nl). We find that the logarithmic eigenvalue repulsion of RMT is exact 
for T^'s close to unity, but overestimates the repulsion for weakly transmit- 
ting channels. The non- logarithmic repulsion resolves several long-standing 
discrepancies between RMT and microscopic theory, notably in the magni- 
tude of the universal conductance fluctuations in the metallic regime, and in 
the width of the log-normal conductance distribution in the insulating regime. 

PACS numbers: 72.10.Bg, 05.60.+W, 72.15.Rn, 73.50.Bk 
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I. INTRODUCTION 



A fundamental problem of mesoscopic physics is to find the statistical distribution of the 
scattering matrix in an ensemble of disordered conductors. Once this is known, one can com- 
pute all moments of the conductance, and of any other transport property, at temperatures 
which are sufficiently low that the conductor is fully phase-coherent. Random-matrix theory 
(RMT) addresses this problem on the basis of the assumption that all correlations between 
the transmission eigenvalues are due to the jacobian from matrix to eigenvalue space .^Q The 
transmission eigenvalues Ti,T2, . . .T^r are the eigenvalues of the matrix product tt\ where 
t is the N X N transmission matrix of the conductor. The jacobian is 

J{{K}) = ll\^j-X^\^ (1.1) 

i<j 

where Aj = (1 — Ti)/Ti is the ratio of reflection to transmission probabilities (A > 0, since 
< T < 1), and 13 G {1, 2, 4} is the symmetry index of the ensemble of scattering matrices. 
[In the absence of time-reversal symmetry, one has /3 = 2; In the presence of time- reversal 
symmetry, one has /? = 1 (4) in the presence (absence) of spin-rotation symmetry.] 

If all correlations are due to the jacobian, then the probability distribution 
P(Ai, A2, . . . Aat) of the A's should have the form P oc JJli f{K), or equivalently, 

P({A„}) = Cexp[-P(j2 uiK A,) + E ^(A.))] , (l-2a) 

i<j i 

u(A,,A,) = -ln|A, -A,|, (1.2b) 

with V = — /3~^ln/ and C a normalization constant. Eq. (|1.2|) has the form of a Gibbs 
distribution at temperature f3~^ for a fictitious system of classical particles moving in one 
dimension in an external potential V, with a logarithmically repulsive interaction u. All 
microscopic parameters (sample length L, width W, mean free path /, Fermi wave length 
Ap) are contained in the single function ^(A). The logarithmic repulsion is independent of 
microscopic parameters, because of its geometric origin. 

The RMT probability distribution (|1.2| ), due to Muttalib, Pichard, and Stone, was jus- 
tified by a maximum-entropy principle for quasi-one-dimensional (quasi-lD) conductors.!'! 
Quasi-ID means L ^ W. In this limit one can assume that the distribution of scatter- 
ing matrices is only a function of the transmission eigenvalues (isotropy assumption). The 
distribution ( p,.2| ) then maximizes the information entropy subject to the constraint of a 
given density of eigenvalues. The function V{X) is determined by this constraint and is not 
specified by RMT. 

It was initially believed that Eq. ( |1.2|) would provide an exact description in the quasi- ID 
limit, if only V^(A) were suitably chosen.ifl However, it was shown recently by one of uM 
that RMT is not exact, even in the quasi- ID limit. If one computes from Eq. (|1.2|) in the 
metallic regime the variance VarG of the conductance G = GoJ^nTn (with Gq = 2e^//i), 
one findU 

VarG/Go = (1.3) 
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independent of the form of V{X). The diagrammatic perturbation theoryoQ of universal 
conductance fluctuations (UCF) gives instead 

VarG'/Go = ^r' (1-4) 

for a quasi- ID conductor. The difference between the coefficients | and ^ is tiny, but it 
has the fundamental implication that the interaction between the A-variables is not pre- 
cisely logarithmic, or in other words, that there exist correlations between the transmission 
eigenvalues over and above those induced by the jacobian. 

What then is the status of the random-matrix theory of quantum transport? It is obvi- 
ously highly accurate, so that the true eigenvalue interaction should be close to logarithmic. 
Is there perhaps a cutoff for large separation of the A's? Or is the true interaction a many- 
body interaction, which can not be reduced to the sum of pairwise interactions? That is 
the problem addressed in this paper. A brief account of our results was reported in a recent 
Letter.! 

The transport problem considered here has a counterpart in equilibrium. The Wigner- 
Dyson RMT of the statistics of the eigenvalues {E^} of a random hamiltonian yields a 
probability distribution of the form ( p.. 21 ), with a logarithmic repulsion between the energy 
levels.§ It was shown by Efetov0 and by Al'tshuler and ShklovskiiLl that the logarithmic level 
repulsion in a small disordered particle (diameter L, diffusion constant D) holds for energy 
separations small compared to the Thouless energy = TiD / . For larger separations the 
interaction potential decays algebraically.0 As we will see, the way in which the RMT of 
quantum transport breaks down is quite different: The interaction u(Aj, A-,) = — In |Aj — Aj| 
is exact for Aj, \j <^ 1, i.e. for strongly transmitting scattering channels (recall that A <C 1 
implies T = (1 + A)~^ close to unity). For weakly transmitting channels the repulsion is 
still logarithmic, but reduced by a factor of two from what one would expect from the 
jacobian. This modified interaction explains the | — ^ discrepancy in the UCF in the 
metallic regime,! and it also explains a missing factor of two in the width of the log- normal 
distribution of the conductance in the insulating regime.lli 

Our analysis is based on the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation 

'|^-^5^^^|:|-A.(l + A.,.Aj-'P. (1.5) 

with ballistic initial condition lim^^o-P = Dj ^(-^j ~ 0+), which describes the evolution of 
the eigenvalue distribution function in an ensemble of disordered wires of increasing length. 
Eq. ( |1.5| ) was derived by Dorokhov,0 (for (3 = 2) and by Mello, Pereyra, and Kumar,El (for 
(3=1, with generalizations to /5 = 2, 4 in Refs. [T^JT^) by computing the incremental change 
of the transmission eigenvalues upon attachment of a thin slice to the wire. It is assumed 
that the conductor is weakly disordered, / 3> Ap, so that the scattering in the thin slice can 
be treated by perturbation theory. A key simplification is the isotropy assumption that the 
fiux incident in one scattering channel is, on average, equally distributed among all outgoing 
channels. This assumption restricts the applicability of the DMPK equation to the quasi- ID 
regime L ^ W, since it ignores the finite time scale for transverse diffusion. 
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Eq. (p..5|) has the form of a diffusion equation in a comphcated A^-dimensional space 



(identified as a certain Riemannian manifold in Ref. 18). For a coordinate-free "supersym- 
metry formulation" of this diffusion process, see Refs. p!9|j20| . The similarity to diffusion in 
real space has been given further substance by the demVnWration§ that Eq. ( |1.5|) holds on 
length scales 3> / regardless of the microscopic scattering properties of the conductor. 

The diffusion equation ( |1.5|) has been studied extensively for more than ten years. Exact 
solutions have been obtained by Mel 'nikovH and Melloil for the case = 1 of a single de- 
gree of freedom (when J = 1). For > 1 the strong coupling of the scattering channels by 
the jacobian ( |1 . 1| ) prevented an exact solution by standard methods. The problem simplifies 
drastically deep in the localized regime {L ^ Nl), when the scattering channels become 
effectively decoupled. Pichardlll has computed from Eq. ( |1.5|) the log-normal distribution 
of the conductance in this regime, and has found an excellent agreement with numerical 
simulations of a quasi-lD Anderson insulator. In the metallic regime {L -C Nl), Mello and 
Stone0@ were able to compute the first two moments of the conductance, in precise agree- 
ment with the diagrammatic perturbation theory of weak localization and UCF [Eq. (p..4[Ll 
in the quasi- ID limit. (Their method of moments has also been applied to the shot noise,E^ 
where there is no diagrammatic theory to compare with.) More general calculations of the 
weak localization effect^ and of universal fluct national [for arbitrary transport properties 
of the form A = J^n^^iTn)] were recently developed, based on linearization of Eq. ( |1.5| ) in 
the fluctuations of the A's around their mean positions (valid in the large- A^ metallic regime, 
when the fluctuations are small). The work of Chalker and MacedocH was motivated by the 
same | — ^ discrepancyi as the present paper and Ref. H, with which it has some overlap. 

None of these calculations suffices to determine the form of the eigenvalue interaction, 
which requires knowledge of the complete distribution function. Here we wish to present (in 
considerable more detail than in our LetteiS) the exact solution of Eq. ( |1.5| ) for /3 = 2. 

The outline of this paper is as follows. In Sec. II we solve Eq. ( p..5|) exactly, for all A^ and 
L, for the case /3 = 2. The method of solution is a mapping onto a model of non-interacting 
fermions, inspired by Sutherland's mapping of a different diffusion equation.!! The case 
P = 2 is special, because for other values of P the mapping introduces interactions between 
the fermions. The free-fermion problem, which is obtained for (3 = 2, has the character 
of a one-dimensional scattering problem in imaginary time. The absence of a ground state 
is a significant complication, compared with Sutherland's problem.c^Ell The exact solution 
which we obtain has the form of a determinant of an A^ x A^ matrix. The determinant can be 
evaluated in closed form in the metallic regime L <^ Nl and in the insulating regime L 3> Nl. 
These two opposite regimes are discussed separately in Sees. Ill and IV. We conclude in Sec. 
V with a comparison of the solution of Eq. ( |1.5| ) with the probability distribution ( |1.2D of 
random-matrix theory. 



II. EXACT SOLUTION 

The solution of the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation (|1.5| ) proceeds in 
a series of steps, which we describe in separate subsections. 
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A. Transformation of variables 



The DMPK equation ( |1.5[ ) can be written in the form of an A^-dimensional Fokker-Planck 
equation, 

D(A) = -A(l + A), (2.1b) 

7 

n{{Xn}) = -J2l^\Xj~Xi\, (2.1c) 

i<j 

where we have abbreviated s = L/l, •y = [3N + 2-/3. Eq. ( p^.lj ) is the diffusion equation in 
"time" s of a one- dimensional gas of classical particles with a logarithmically repulsive 
interaction potential Vt. The diffusion takes place at temperature (3~^ in a fictitious non- 
uniform viscous fluid with diffusion coefficient D{X). 

The position dependence of the diffusion coefficient is problematic. We seek to 
eliminate it by a transformation of variables. Let {x„} be a new set of N indepen- 
dent variables, related to the A's by A„ = f{xn)- The new probability distribution 
P{{xn},s) = P{{Xn}, s)Yli\f'{xi)\ stiU satisfies a Fokker-Planck equation, but with a 
new potential fl{{xn}) and a new diffusion coefficient D{x). The potential transforms as 
f2 — i> f2 — /3~^I]jln \f'{xi)\, while the diffusion coefficient transforms as D —>■ D/ f'{xY. In 
order to obtain an x-independent diffusion coefficient, we thus need to choose f{x) such that 
/(x)[l + /(a;)]//'(a;)^ = constant. The choice f{x) = sinh^ x does it. 

We therefore transform to a new set of variables defined by 



A„ = sinh Xn, T„ = 1/ cosh x„. (2.2) 

Since T„ G [0, 1], x„ > 0. The probability distribution of the x- variables satisfies a Fokker- 
Planck equation with constant diffusion coefficient, 

,) ^ i g A + ffpirn(M)) , (2.3a) 

Q{{xn}) = — ^ In I sinh^ Xj — sinh^ Xj| — — ^ In | sinh2xj|. (2.3b) 



i<j 



P 



It turns out that the x- variables have a special physical significance: The ratio L/x„ equals 
the channel-dependent localization length of the conductor.! 



B. Prom Fokker-Planck to Schrodinger equation 

Sutherlandii has shown that a Fokker-Planck equation with constant diffusion coefficient 
and with a logarithmic interaction potential can be mapped onto a Schrodinger equation 
with an inverse-square interaction which vanishes for (3 = 2. The Fokker-Planck equation 
^ does have a constant diffusion coefficient, but the interaction is not logarithmic. It 
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is not obvious that Sutherland's mapping onto a free-fermion problem should work for the 
non-translationally invariant interaction (p.3b| ), but surprisingly enough it does. 

To map the Fokker-Planck equation ( |2.3| ) onto a Schrodinger equation we substitute 

P({xJ, s) = exp [-lf3n{{xn})] ^{{xn}, s) . (2.4) 

This is a variation on Sutherland's transformation,ii which we used in Ref. ^ in a different 
context. Substitution of Eq. ( |2.4| ) into Eq. ( |2.3a| ) yields for the equation 



9^ 1 ^d^^ P 
— = > h— \E'> 

ds 27^9x2 47 



2 \ dxi I 



dxi 



(2.5) 



The expression between square brackets is evaluated as follows (we abbreviate = cosh 2xi): 



4EE 



4 



'2 



- 6 



1 



4 



P - \ 4 

p -1 

4E E E — - — 

P — 1 4 1 



(2.6) 



+ -E72^ 



In the final equality we have used that for any three distinct indices i,j, k 



. 3 ; 



(2.7) 



+ 



+ 



i^j - 6) i^k - e.) (6 - Qi^k - Q (e. - ^k) (e, - e 



1, 



(2.8) 



so that the triple sum over k ^ i ^ j collapses to a double sum over i ^ j. Collecting 
results, we find that satisfies a Schrodinger equation in imaginary time, 



(9^ 

ds 

n 



{n-u)m, 



27 



sinh 2xi 



I /?(/?-2) ^ 



27 



sinh^ 2a;j + sinh^ 2xi 
(cosh2a;j — cosh2a;i)" 



U = N(N - 1)- - N(N - 1)(N - 2)^. 

27 7 67 



(2.9a) 
(2.9b) 

(2.9c) 



The interaction potential in the hamiltonian (|2.9bD is attractive for /3 = 1 and repulsive 
for /3 = 4. For /? = 2 the interaction vanishes identically, reducing 7i to a sum of single- 
particle hamiltonians TYq, 



Hp, 



1 



1 



ANdx"^ 4A^sinh2 2a;' 
6 



(2.10) 



(Note that 7 = 2N for (3 = 2.) It might be possible to solve also the interacting Schro- 
dinger equation ( p.9|) ioiB = 1 or 4, by some modification of techniques developed for the 
Sutherland hamiltonianjH§Ey but in this paper we focus on the simplest case /5 = 2 of broken 
time-reversal symmetry. 

To complete the mapping onto a single-particle problem, we need to consider the bound- 
ary condition at the edge a; = 0. (Recall that x > 0.) Conservation of probability implies 
for P the boundary condition (one for each i = 1, 2, . . . A^) 

>-fl^+ ^^1^1-0. (2.11) 



^0 Y dxi dxi J 

According to Eq. (p.4|), the corresponding boundary condition on ^' is 



+0 \dxi ^ dx 
which in view of Eq. (|2.3b[ ) simplifies to 



lim (1^ - . ^ ] = 0, (2.13) 



*o \ dxi sinh 2x,- 



independent of /?. Fortunately, the boundary condition does not couple different degrees of 
freedom, so that we have indeed obtained a single-particle problem for (3 = 2. 



C. Prom probability distribution to fermion Green's function 

We seek a solution P({x„}, s | of the Fokker-Planck equation (|2.3| ) with sym- 

metrized delta-function initial condition 

1 ^ 

P{{xn], I {vn}) = Try E n ^(^^ - (2-14) 

■ TT i = l 



The sum in Eq. ( p.l4| ) is over all A^! permutations of 1, 2, ... iV. Eventually, we will take the 



limit {Un} ^ of a ballistic initial condition, but it is convenient to first consider the more 
general initial condition (|2.14|). In this subsection we use the mapping onto a Schrodinger 



equation of the previous subsection to relate the probability distribution P({x„},s | {Vn}) 
to the iV-fermion Green's function G{{xn\,s \ {Un})- 

We first note that, since exp(— is an s- independent solution of the Fokker-Planck 
equation (^73|) , exp(— is an s-independent solution of the Schrodinger equation ( p.9|) 
[in view of the mapping ( p.4| )]. For a particular ordering of the x„'s, the function \l/o oc 
exp(— is therefore an eigenfunction of the A^-fermion hamiltonian Ti with eigenvalue 
U . Anti-symmetrization yields the fermion eigenstate 

^o({x4) = C^exp \-\(3n{{xM n ^l—^ (2.15) 



with C a normalization constant. 
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We obtain the A^-fermion Green's function G from the probabihty distribution P by the 
similarity transformation 



G{{Xn}. S I {yn}) = %\{Xn})Pi{Xn}, S \ {Vn})^ .{{Vn}) ■ (2.16) 

To verify this, we first observe that G is by construction anti-symmetric under a permutation 
of two X or two y variables. For a given order of the XnS, the function G satisfies the 
Schrodinger equation 

-^ = {n-U)G, (2.17) 



in view of Eqs. ( |2.4| ), ( p.9a[ ), and ( |2.15| ). Finally, Eq. ( |2.14| ) implies the initial condition 



1 ^ 

G{{xn}, I {y^}) = ^ E n -y^.)^ (2-18) 

■ i = l 

with (Ttt the sign of the permutation. Hence G is indeed the iV-fermion Green's function. 
The relation ( p^.l6 ) holds for any j3. In the remainder of this paper we consider the 



non-interacting case (3 = 2. The eigenstate (|2.15|) then takes the form 

^^{{xn}) = Cn(sinh^a;j -sinh2xi)n(sinh2xi)^/^ (2.19) 

i<j i 

The iV-fermion Green's function G becomes a Slater determinant of the single-particle 
Green's function Go, 

G{{xn}, s I {yn}) = -j^ Det Go{xn, s\ym), (2.20) 

where Det a^.^ denotes the determinant of the NxN matrix with elements a„m- The function 
Go{x, s\y) is a. solution of the single-particle Schrodinger equation —dGo/ds = TCqGq in the 
variable x, with initial condition G{x, 0\y) = 6{x — y)- In the following subsection we will 
compute the single-particle Green's function Gq. The probability distribution P, for (3 = 2, 
then follows from Eqs. (pl6| ), (^J^), and ( p:20| ): 



^/r ,r IN n»vi(sinh^x,- — sinh^Xj) n»(sinh2a:i)-'^/^ e^* _ ^ , , . , 

^«^»'- ^ I (^"» = lCW^-^mh-Jn'U2j-/. 1^ Go(.„. . I (2.21) 



D. Computation of Green's function 

To compute the Green's function Gq of the single-particle hamiltonian ( |2.10| ) we need to 
solve the eigenvalue equation 

1 ,, , 1 iIj(x) ,, , , , 



with the boundary condition dictated by Eq. ( ^.13| ) 
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We have found that the substitution 

^(x) = (sinh2x)^/V(cosh2a;) (2.24) 
transforms Eq. ( |2.22| ) into Legendre's differential equation in the variable z = cosh2x, 

-.{Ne + l)fiz). (2.25) 



d_ 

dz 



The boundary condition ( |2.23| ) restricts the solutions of Eq. ( p.25| ) to the Legendre functions 
of the first kind Vy[z). The index u is given by z/ = — ^ + |iA; with k a real number. (These 
Legendre functions are also known as "toroidal functions" , because they appear as solutions 
to the Laplace equation in toroidal coordinates.) The numbers k, and e are related by 
— z/(z/ + 1) = Ne + i and e = ^k'^/N. We can restrict ourselves to > 0, since the functions 
P 1 , iji. and P are identical. 

We conclude that the spectrum of TCq is continuous, with positive eigenvalues e = \k'^ /N. 
The eigenfunctions ipk{x) are real functions given by 

ipkix) = [7rA;tanh(i7rfc)sinh(2x)]^/^Pi(ifc_i)(cosh2x). (2.26) 

They form a complete and orthonormal set, 

dk tlJk{x)Mx') = 27r(5(a; - x'), (2.27) 

oo 

dxtlJk{x)^k'{x) = 2n6{k - k'), (2.28) 











in accordance with the inversion formula in Ref. The single-particle Green's function Gq 
has the corresponding spectral representation 



Go{x, s\y) = (27r)-i / dk expi-l^s/N) Mx)My) 
Jo 

= |(sinh2a;sinh2?/)^/2 / dk exp{-\k^s/N)ktanh{^TTk) 

X Pi(i,_i)(cosh2x)P.(i,_i)(cosh2y). (2.29) 

E. Ballistic initial condition 

Eqs. ( |2.21| ) and ( f^.29D together determine the probability distribution P({x„},s | {i/n}) 
with initial condition (|2.14| ), 
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P = Cis[ 



ni<j(sinh^ Xj — sinh^ Xi) nj(sinh 2xj 



ni<j(sinh yj - sinh yi 



X Det 



roo 

/ dk exp(— i/c^s/A^)A;tanh(|7rA;)Pi(j^_;,^-)(cosh2x„)Pi(i^_;^)(cosh2?/m) 
C(s) ]^(sinh^ Xj — sinh^ Xj) PJ(sinh 2xi) j dki J dk2 ■ ■ ■ J dk^ 



Y[ \_QW{—jkfs/N)ki tanh(i7rA;j)Pi(;^._^-)(cosh 2x 



Det Pi(ifc^_i) (cosh 2ym] 
Ui<j (sinh'^ yj - sinh^?/i 



(2.30) 



We have absorbed all x and y independent factors into the function C{s), which is fixed by 
the requirement that P is normalized to unity, 



roo roo roo 

/ dxi / dx2 ■ ■ ■ / dxN P = 1. 

Jo Jo Jo 



(2.31) 



In the second equality in Eq. ( p.30|) we have applied the identity Det {hnanm) = 
(Jli bi)Detanm to isolate the factors containing the y-variables. 

Now it remains to take the limit {y-n} — of a ballistic initial condition. The limit is 
tricky because it involves a cancellation of zeroes of the determinant in the numerator with 
zeroes of the alternating function in the denominator. It is convenient to first write the 
alternating function Vandermonde determinant, 



JKsinh^ yj - sinh^ yi) = Det (sinh^ y^T ^■ 

i<j 

Next, we expand the Legendre function in powers of sinh^y, 

oo 

p.(i,.i)(cosh2y) = Y.c,m^m'^'yY-'- 
p=i 

The factors Cp{k) are polynomials in fc^, with Ci{k) = 1 and 

Cpik) = i-ir-'[2P~\p - l)!]^2 (^2 ^ ^2)(^2 + 32) . . . ^^2 ^ - 3)=^) 



(2.32) 



(2.33) 



(2.34) 



for p > 2. In the limit y — > 0, we can truncate the expansion ( p. 331) after the first N terms, 
that is to say. 



lim 



DetP 



i(ifc„-l) 



(cosh 2?/„ 



lim 



Det 



E^=iCp(A;„)(sinh2y^)f-i 



{?/n}^o J].^j.(sinh yj — sinh y^) {yn}^o 



Det (sinh y^ 



\n-l 



(2.35) 



The numerator on the r.h.s. of Eq. (|2.35| ) factors as the product of two determinants, one of 
which is just the Vandermonde determinant in the denominator, so that the whole quotient 
reduces to the single determinant DetCm(fcn)- This determinant can be simplified by means 
of the identity 



Det Cmikn) = Co Det (fc, 



2\m-l 



(2.36) 
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with Co a numerical coefficient. Eq. ( |2.36|) holds because the determinant of a matrix is 
unchanged if any one column of the matrix is added to any other column, so that we can 
reduce the polynomial Cm{k) in k"^ of degree m — 1 to just its highest order term /^^(m-i) 
times a numerical coefficient. 
Collecting results, we find 



DetPi(i;,„_l)(cOSh2yJ ^w,2^m-l 

hm — — ^\ = coDet (A;„)"' \ 



(2.37) 



Substituting into Eq. ( p. 301) , and absorbing the coefficient Cq in the function C(s), we obtain 
the probability distribution P{{xn}, s) for a ballistic initial condition. 



P{{xn}, s) = C{s) ]^(sinh^Xj — sinh^Xj) ]^(sinh2a;i) 



i<j 



X Det 



dk exp{-lk'^s/N) tanh(|7rfc)A;2'"-^ (cosh 2x„ 



(2.38) 



This is the exact solution of the DMPK equation for the case j3 = 2. 



III. METALLIC REGIME 
A. Probability distribution 



The solution (|2.38| ) holds for any s and A^. It can be simplified in the regime 1 ^ s <^ 



of a conductor which is long compared to the mean free path I but short compared to 
the localization length Nl. This is the metallic regime. The dominant contribution to the 
integral over k in Eq. ( |2.38D then comes from the range k ^ {N/sY^"^ ^1. In this range 



tanh(i7r/c) 1 and the Legendre function simplifies to a Bessel function,! 

Pi(ifc_i)(cosh2x) = Joikx) ( .^^ ^^ j for A; > 1, x > 1/k. (3.1) 

The second condition x ^ 1/k on Eq. (|3.1|) implies the restriction x ^ (s/A^)^/^, which is 
irrelevant since s/N ^ 1. The fc-integration can now be carried out analytically. 



/■oo 

/ ^ 

Jo 



dk expi-^kh/mk^""-^ Mkxn) = ^(m - 1)! (4N/sr exp(-xlN/s)Lm-i(xlN/s), 



(3.2) 

with Lm-i a Laguerre polynomial. We then apply the determinantal identity 

BetLm-i{xlN/s) = cDet (xl^'^ = cJlla^l - x^), (3.3) 

i<j 

with c an x-independent number [which can be absorbed in C{s)]. Eq. ( p.3| ) is derived in the 
same way as Eq. ( p.36|) , by combining columns of the matrix of polynomials in x^. Collecting 
results, we find that the general solution ( ^2.38] ) simplifies in the metallic regime to 
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P{{xn},s) = C{s)Y\_ [(sinh^Xj — sinh^Xj)(x^ ~ 2;^)] H [6xp(— x^A^/s)(xj sinh 2xj)"^ 



/2' 



(3.4) 



In the remainder of this section we use the probabihty distribution ( |3.4| ) to compute 
various statistical quantities of interest. For that purpose it is convenient to write P as a 
Gibbs distribution, 

P{{xn}, s) = C{s) exp -l3 (J2 "^i^h Xj) + J2 ^(^»' ' (3-5a) 

i<j i 

u{xi, Xj) = — I In I sinh^ Xj — sinh^ Xi| — | In |xj — x^l, (3.5b) 

V{x, s) = ^Ns-^x^ - i ln(x sinh 2x), (3.5c) 

with /3 = 2 understood. 



B. Eigenvalue density 

The mean density {p{x))s of the a;-variables is defined as the ensemble average with 
distribution P({x„},s) of the microscopic density p{x): 

N 



n=l 

{p{x))s= dxi dx2---j dxN P{{xn},s)p{x). (3.7) 
Jo Jo Jo 

The mean density is determined to leading order in by the integral equation 

— / dx' {p{x'))su{x, x') = V{x, s) + const. (3.8) 



The additive constant (which may depend on s but is independent of x) is fixed by the 
normalization condition 

/•oo 

dx {p{x))s = N. (3.9) 

Eq. ( ^.81 ) can be understood intuitively as the condition for mechanical equilibrium of a fic- 
titious one-dimensional gas with two-body interaction m in a confining potential V. Dysonil 
has shown that corrections to Eq. ( p.8[ ) are an order A^~^ In smaller than the terms re- 
tained, and are /3-dependent. These corrections are responsible for the weak-localization 
effect in the conductance.^ Here we consider only the leading order contribution to the 
density, which is of order N and which is independent of (3. 

Substituting the functions u{x, x') and V{x, s) from Eq. ( |3.5| ) into Eq. ( p.8| ), and taking 
the derivative with respect to x to eliminate the additive constant, we obtain the equation 

s /""^r / , / / sinh 2a; 2x \ ,^r-. 

— / dx' {p{x'))s —2^ + ^ ]=x + 0{1N). 3.10 

2i\ Jo \ smh X — smh x' x^ — x''^ 
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We note that 



dx' 



sinh 2x 



sinh X — sinh x' 



+ 



2x 



x" 



X' 



12 



In 



sinh(s + x) 



sinh(s 



X] 



+ In 



s + X 



X 



2x + 0{x/s), for s > 1, s > x. 



(3.11) 



It follows that the uniform density 



(P(^)). 



N 



9{s~x) 



(3.12) 



is the solution of Eq. ( p.lO| ) in the regime s ^ 1, s ^ x. (The function 6{^) equals 1 for 
^ > and for ^ < 0.) The result ( 3.12| ) was first obtained by Mello and Pichard, by direct 
integration of the DMPK equation.^ To order N , the x-variables have a uniform density 
of Nl/L, with a cutoff at L/Z such that the normalization ( p. 91) is satisfied. In the cutoff 
region x ^ L/l the density deviates from uniformity, but this region is irrelevant since the 
transmission eigenvalues are exponentially small for x ^ 1. 



C. Correlation function 

The two-point correlation function K{x,x', s) is defined by 

K{x, x', s) = {pix))s {pix'))s - {p{x)p{x'))s. (3.13) 

We compute the two-point correlation function by the general method of Ref. which 
is based on an exact relationship between K and the functional derivative of the mean 
eigenvalue density (p) with respect to the eigenvalue potential V: 



Eq. ( p. 14 ) holds for any probability distribution of the form ( p.5a ), regardless of whether 



the interaction is logarithmic or not. In the large- iV limit the functional derivative can be 
evaluated from the integral equation (|3.8|). The functional derivative S{p)/SV equals the 
solving kernel of 



— / dx' tp{x')u{x, x') = (f){x) + const, (3.15) 
Jo 

where the additive constant has to be chosen such that ip has zero mean, 

dx'^{x) = 0, (3.16) 



since the variations in (p) have to occur at constant A^. Because of Eq. ( p. 14 ), the integral 
solution 

i){x)= dx' pK{x,x')4>{x') (3.17) 
Jo 
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of Eq. ( |3.15|) directly determines the two-point correlation function. It turns out that 
K{x, x') = K{x, x', s) is independent of s in the metallic regime. 

The integral equation ( p.l5| ) can be solved analytically by the following method. We 
extend the functions ip and symmetrically to negative x, by defining ■?/'(— x) = ipix), 
4>{—x) = 4>{x). We then note that the decomposition 

u{x, x') = -\U{x - x') - \ U{x + x) + In 2, (3.18a) 
= In |2a; sinh x|, (3.18b) 



transforms the integral equation (|3.15| ) into a convolution, 

D 

dx' i){x')U{x - x') = 2(f){x) + const, (3.19) 
which is readily solved by Fourier transformation. The Fourier transformed kernel is 

dxe'^^'Uix) = ---[1 + cotanh(i7r|fc|)]. (3.20) 

-oo I rC I 

The fc-space solution to Eq. ( |3.19|) is ipi^k) = 2(j){k)/U{k), which automatically satisfies the 
normalization ( p.l6|) . In x-space the solution becomes 

/oo 
dx' }C{x - x')(j){x') 
-oo 
rca 

= 2/ dx'\]C{x-x') + )C{x + x')](f){x'), (3.21a) 
Jo 

/C(x) = — / dke = - dk——. (3.21b) 

^ ' 27ri-oo U{k) n Jo U{k) ^ ' 

Combining Eqs. ( p.l7|) , (|3.2CI|) , and (|3.21| ), we find that the two-point correlation function 
is given by 

K{x,x') = }C{x-x') + IC{x + x'), (3.22a) 

JC{ ) ^ [°°dk ^cosfcx (3 22b) 

Pn"^ Jo 1 + cotanh(|7rA;) ' 

with (3 = 2. The inverse Fourier transform ( p.22b|) evaluates to 

1 



^ Re 



fx + iO+)"^ - (x + iTT)"^ 



(3.22c) 



where O"*" is a positive infinitesimal. 

We derived i these expressions for the two-point correlation function for the case (3 = 2. A 
direct integration of the DMPK equation by Chalker and Macedd23 shows that the function 
K{x,x') has in fact the 1//? dependence indicated in Eq. (|3.22|) , as expected from general 
considerations .0 
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D. Universal conductance fluctuations 



Now that we have the two-point correlation function, we can compute the variance 
Var A = (A^) — (A)'^ of any hnear statistic A = J2n=i o-i^n) on the transmission eigenvalues 
(recall that T„ = cosh~^x„). By definition 



VarA = -/ dx dx' a{x)a{x')K{x,x'). (3.23) 
JO Jo 

Substituting Eq. (^1221) we find 

VarA = -— / dk ' ^ ;' 3.24a 

I3tt^ Jo 1 + cotanh(i7rA;)' ^ ^ 

a{k) = 2 / dxa{x) cos kx, (3.24b) 



or equivalently, 

Var A = -- Hd. Hd.' ( ( 'JP-] lu ( i±4(ii±4^') . (3.24c) 

2(371^ Jo Jo \ dx J \ dx' J \l + 7T^{x-x')-^J ^ ' 

To obtain the variance of the conductance G/Gq = I]„T„ (with Go = 2e^//;.), we substi- 
tute a{x) = cosh~^a;, hence a{k) = nk/ sinh(i7rfc), hence 

VaTG/Go = (3.25) 
15 

in agreement with Eq. ( [1.41 ). In the same way one can compute the variance of other 
transport properties. For example, for the shot-noise power!! P/Po = Z]n^n(l — ^n) (with 
Pq = 2e\V\Go and V the applied voltage) we substitute a{x) = cosh~^x — cosh'^'x, hence 
a{k) = ^TTk(2 — A;^)(sinh ivrA;)"^, hence 

46 

VarP/Po = ^r\ (3.26) 

in agreement with the result obtained by a moment expansion of the DMPK equation.!! 
Another example is the conductance Gns of a normal-superconductor junction, which for f3 = 
1 is a linear statistic,!! Gns/Gq = En2T2(2 - T„)-2. We substitute a{x) = 2cosh"^a; (2 - 
cosh~^a;)^ = 2 cosh~^(2a;), hence a{k) = ^wk/ sinh(|7r/c), hence 

VarGNs/Go = ^-^. (3.27) 

15 TT* 

Finally, for the variance of the critical current Ic of a point-contact Josephson junction 
(which is also a linear statistic for /? = 1)011 we compute 

Var Jc//o = 0.0890, (3.28) 

with /q = eA/h and A the superconducting energy gap. 

As in the previous subsection, we note that our results are derived for /3 = 2, and that 
the 1/P dependence of the variance in Eq. ( p.24| ) needs the justification provided by the 
calculation of Chalker and Macedo.i! 
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IV. INSULATING REGIME 



The solution (|2.38|) can also be simplified in the regime 1 -C <^ s of a conductor which 
is long compared to the localization length Nl. This is the insulating regime. It is sufficient 
to consider the range x„ ^ 1, since the probability that x ^ 1 is of order N/s which is <^ 1. 
The appropriate asympotic expansion of the Legendre function is 



Pi(ifc-i) (cosh 2x) = {2n sinh 2x)-^^^ Re for x > 1. (4.1) 

For s/N ^ 1, the dominant contribution to the integral over k in Eq. (|2.38|) comes from 
the range -C 1. In this range tanh(i7rfc) — >■ ^7rk and the ratio of Gamma functions in Eq. 
( [4.1|) simplifies to 

= i^k^^y fol- < 1- (4-2) 



r(i + iiA;) 

The /c-integration can now be carried out analytically, 

J^^dk exp{-lkh/N)k^^-^ sm{kxn) = i-ir~^TT'^\N/sr exp{-xlN/s)H2rn-i{xn\/N/~s), 

(4.3) 

with H2m-i a Hermite polynomial. We then apply the determinantal identity [cf. Eq. (p. 3D] 



Deti/2m-i(x„ViV/s) = cDetx^"-^ = cl[x, l[{x] - x^), (4.4) 

i i<j 



with c an x-independent number. Collecting results, we find that the general solution (p. 38 
reduces in the insulating regime to 

P{{xn}, s) = C{s) Yl [(sinh^Xj - sinh^Xi)(xJ - x-)] Y[ [exp(-x- A^/s)xi(sinh2xi)^/^ . 

i<j i 

(4.5) 

(This formula was cited incorrectly in our Letter.!) 

The result ( [4.5|) can be simplified further by ordering the x„'s from small to large and 
using that 1 ^ Xi ^ X2 ^ • ■ ■ ^ x^r in the insulating regime {s ^ N). The distribution 
function then factorizes, 

N 

P({xJ, s) = C{s) n exp [{2i - l)x, - x]N/s 

i=l 

N 

= (7rs/iV)-^/2 nexp [-(iV/s)(x, - x,)^] . (4.6) 

i=l 

The x„'s have a gaussian distribution with mean x„ = \{s/N){2n — 1) and variance \s/N . 
The width of the gaussian is smaller than the mean spacing by a factor (N/sY^"^, which is 
<C 1, so that indeed 1 <c; xi -C X2 <^ ■ ■ ■ -C xat, as anticipated. 
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The conductance G/Gq = J2n cosh.~'^ Xn is dominated by Xi, i.e. by the smallest of the 
XnS. Since Xi ^ 1 we may approximate G/Gq = 4exp(— 2a;i). It follows that the con- 
ductance has a log-normal distribution, with mean (InG/Go) = —s/N + 0{1) and variance 
Var InG/Go = 2s/N. Hence we conclude that 

Var InG/Go = -2(lnG/Go), (4.7) 

in agreement with the result obtained by Pichard,lil by directly solving the DMPK equation 
in the localized regime. 

The results obtained here are for the case (3 = 2. Pichard has shown that the relationship 
( [4.7| ) between mean and variance of InG/Go remains valid for other values of (3, since both 
the mean and the variance have a dependence on the symmetry index. 



V. COMPARISON WITH RANDOM-MATRIX THEORY 

The random-matrix theory of quantum transportiH is based on the postulate that all 
correlations between the transmission eigenvalues are due to the jacobian ( |1.1D . The resulting 
distribution function ( p..2| ) has the form of a Gibbs distribution with a logarithmic repulsive 
interaction in the variables A„ = (1 — T„)/T„. There exists a maximum-entropy argument 
for this distribution,!'! but it has no microscopic justification. In this paper we have shown, 
for the case of a quasi-lD geometry without time-reversal symmetry, that the prediction of 
RMT is highly accurate but not exact. 

In the metallic regime {L <^ Nl), the distribution is given by Eq. (|3.5| ). In terms of 
the A- variables (A = sinh^x), the distribution takes the form (|1.2a|) of RMT, but with a 
different interaction 

u{Xi,Xj) = — |ln|Aj — Aj| — I In larcsinh^A^''^ — arcsinh^Aj^''^|. (5.1) 

For A <^ 1 (i.e. for T close to unity) u{Xi, Xj) — > — In \Xj — Aj|, so we derive the logarithmic 
eigenvalue repulsion ( |1.2b|) for the strongly transmitting scattering channels. However, for 
A ~ 1 the interaction ( p.l| ) is non-logarithmic. For fixed Aj ^ 1, u{Xi, Xj) as a function of 
Xj crosses over from — In |Aj — Ai| to — | In \ Xj — Xi\ as Aj oo (see Fig. 1). It is remarkable 
that, for weakly transmitting channels, the interaction is twice as small as predicted by 
considerations based solely on the jacobian. We have no intuitive argument for this result. 
The reduced level repulsion for weakly transmitting channels enhances the variance of the 
conductance fluctuations above the prediction ( |1.3| ) of RMT. Indeed, as shown in Sec. III.D, 
a calculation along the lines of Ref . |^, but for the non-logarithmic interaction ( ^.1|) , resolves 
the I — ^ discrepancy between RMT and diagrammatic perturbation theory, discussed in 
the Introduction. The discrepancy is so small because only the weakly transmitting channels 
(which contribute little to the conductance) are affected by the non-logarithmic interaction. 

In the insulating regime {L ^ Nl), the distribution is given by Eq. ( [4.5|) . In terms of 
the A's the distribution takes the form ( [L.2aD of RMT, but again with the non-logarithmic 



interaction ( |5.1|) . Since In A S> 1 in the insulating regime, the interaction ( |5.1|) may be 
effectively simplified to u{Xi, Xj) = — | In \Xj — Aj|, which is a factor of two smaller than the 
interaction ( |1.2b| ) predicted by RMT. This explains the factor-of-two discrepancy between 
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the results of RMT and of numerical simulations for the width of the log-normal distribution 
of the conductance:Ei RMT predicts Var IyiG/Gq = — (InG/Go)? which is twice as small as 
the correct result ( |4.7| ). 

We conclude by mentioning some directions for future research. We have only solved 
the case /3 = 2 of broken time-reversal symmetry. In that case the DMPK equation ( [L.5| ) 
can be mapped onto a free-fermion problem. For /3 = 1,4 the Sutherland-type mapping 
which we have considered is onto an interacting Schrodinger equation. It might be possible 
to solve this equation exactly too, using techniques developed recently for the Sutherland 
hamiltonian.E3E3 From the work of Chalker and Macedoil we know that the two-point cor- 
relation function in the large- limit has a simple l/jS dependence on the symmetry index. 
This poses strong restrictions on a possible /3-dependence of the eigenvalue interaction. 

Another technical challenge is to compute the exact two-point correlation function 
K{x, x', s) from the distribution function P({x„}, s). Our result ( p.5|) for P is exact, but the 
large- asymptotic result ( p.22| ) for K ignores fine structure on the scale of the eigenvalue 
spacing. (This large- iV result for K corresponds to the regime of validity of the diagram- 
matic perturbation theory of UCF,i'0 while the exact result for P goes beyond perturbation 
theory.) In RMT there exists a technique known as the method of orthogonal polynomials JhI 
which permits an exact computation of A logarithmic interaction seems essential for 
this method to work, and we see no obvious way to generalize it to the non-logarithmic 
interaction (|5.1| ). 

It might be possible to come up with another maximum-entropy principle, different from 
that of Muttalib, Pichard, and Stone,i which yields the correct eigjenvalue interaction ( ^.1|) 
instead of the logarithmic interaction ( |1.2b| ). Slevin and Nagao^y have recently proposed 
an alternative maximum-entrooy principle, but their distribution function does not improve 
the agreement with Eq. ( |1.4D .Ej 

To go beyond quasi-one-dimensional geometries (long and narrow wires) remains an 
outstanding problem. A numerical study of Slevin, Pichard, and Muttalib@ has indicated a 
significant break-down of the logarithmic repulsion for two- and three-dimensional geometries 
(squares and cubes). A generalization of the DMPK equation ( |1.5|) to higher dimensions 
has been the subject of some recent invest igat ions. It remains to be seen whether the 
method developed here for Eq. ( |1.5|) is of use for that problem. 
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